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Abstract. Experimental particle spectra can be successfully described by power-law tailed energy distribu- 
tions characteristic to canonical equilibrium distributions associated to Renyi's or Tsallis' entropy formula 
- over a wide range of energies, colliding system sizes, and produced hadron sorts. In order to derive its 
evolution one needs a corresponding dynamical description of the system which results in such final state 
observables. The equations of relativistic fluid dynamics are obtained from a non-extensive Boltzmann 
equation consistent with Tsallis' non-extensive q-entropy formula. The transport coefficients like shear 
viscosity, bulk viscosity, and heat conductivity are evaluate based on a linearized collision integral. 
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1 Introduction 

Experimental particle spectra at the Relativistic Heavy- 
Ion Collider (RHIC) and at the Large Hadron Collider 
(LHC) can be successfully described by power-law tailed 
energy distributions characteristic to canonical equilib- 
rium distributions associated to Renyi's or Tsallis' entropy 
formula - over a wide range of energies, colliding system 
sizes and produced hadron sorts [Tl[^|^l|il[Hl[r)l[71|Rl[^pi71P1 
fT2lfT3liT4lfT5] . Such a generalized statistical model can be 
solicited on the valence quark level, since several experi- 
mental findings on the azimuthal flow component V2 [16L 
[TTlfTS] as well as on the non-extensivity measure (q — 1) 
show quark number scaling [9l [TTfllTl"l[l"2"] . Therefore it is of 
primary importance to establish a theoretical tool for the 
description of the relativistic flow of such a quark matter. 

The non-cxtensivity effects in general can be viewed as 
an effective measure of finite available phase-space |19ll20] , 
either due to spatial finiteness of the reaction zone or due 
to long range entanglement of quark and gluon strings in 
the hadron formation process, representing a leading order 
deviation between microcanonical and canonical approxi- 
mations in traditional terms [20,21 . Taken a hadronizing 
quark matter with effective, massive and entangled quarks 
seriously, the (q— 1) measure for mesons should be about 
half, for baryons and antibaryons about one third of the 
quark value . It means that if the power-law tailed ob- 
served pT-distribution of hadrons was preformed practi- 
cally in the quark phase, then the non-extensivity effects 
should have been even more pronounced on the quark than 
on the hadron level. Therefore it is unavoidable to test 



those models, which are based on or at least are consis- 
tent with non-extensive thermodynamics, deriving their 
equilibrium state based on the Tsallis' or Renyi entropy. 

An important step towards this goal is the develop- 
ment of a relativistic fluid dynamics, which is fully con- 
sistent with non-extensivity. Although initial steps were 
taken by considering the so-called relativistic perfect q- 
hydrodynamics 22 , 23 , 24 , the relativistic equations of mo- 
tion and transport coefficients for a finite q have only re- 
cently been derived [25] . 

In this paper we present the first full derivation of a 
consequent fluid dynamics from a non-extensive relativis- 
tic Boltzmann equation (NEBE). This equation retains 
the form but generalizes the well known Boltzmann trans- 
port equation, with a single parameter denoted by q. We 
derive the relativistic fluid-dynamical equations of motion 
from this non-extensive relativistic Boltzmann transport 
equation. Applying well known traditional methods, we 
calculate the transport coefficients like heat-conductivity, 
shear and bulk viscosity. We also analyse the question 
which entropy and conserved Noether four-currents have 
to be used in a proper, consistent non-extensive approach. 
Finally we also discuss Grad's method of moments to de- 
rive relativistic causal fluid dynamics in this framework. 



2 Non-extensive statistics, kinetic and fluid 
dynamical equations 

A non-extensive ^-generalization of the Boltzmann-Gibbs 
(BG) entropy was proposed by C. Tsallis [TW2rj] based on 
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the following formula, 



w , 

S q = -k B ^2 ln ? (P») = T 
i=i 



,(1) 

where fc^ is the Boltzmann constant, Pi is the probability 
of the ith from a total of W possible microstates, such that 
Pi 6 [0, 1] with the normalization condition ^2%=iPi = 1- 
Here the g-parameter is a real number q g [0, 2], and ln g 
is dehned as the ^-logarithmic function, 



ln 9 (x) = 



A-q _ i 
1-9 



(2) 



for any x > 0. The inverse of the g-logarithmic function is 
the q-exponential, 

exp q (x)^e q (x) = [l + (l-q)x] 1/{1 - q) , (3) 

for any x > 1/ (q— 1). These functions return the well 
known natural logarithm and exponential functions for 
q = 1, that is ln g= i (x) — In (a;) and exp ?=1 (x) = exp (x). 

For later reference we also define the following dual 
functions [2"7] , 



In 



p 9-l 



9* 0) 
9* 0) 



9-1 
[1 + (? - 



l)x 



1/(9-1) 



and 



Ln 9 (x) = q 
E q (x) = 



1 



c^ 1 - 1 

9-1 

g-1 
9 



1/(9-1) 



(4) 
(5) 

(6) 
(7) 

2 — q that 



One can show using Eqs. (I2I3[) with q* 
hi g . (x) = ln 2 _q (x), e r (x) = e 2 - q (x). Similarly, E q (x) = 
e q * (x/q), e q (x) = E 9 » (q*x) and Ln,j (x) = -qln q (1/x) = 
ging* (x). These relations are called dualities and using 
these functions we can re-write the q-generalized entropy 
in two further equivalent forms: 



Ti- 



ll' 



S q = -k B ^piln q , {p^ 



^2p l Ln q (pi 



(8) 



i=i 



We note that the Renyi formula, Sr — — In J2i=i Pi 1S 
preferred in order to simply fulfill the additivity property 
inherent in the zeroth law of thermodynamics [35]. Since 
the Renyi entropy is proportional to the logarithm of a 
phase-space integral, its local four-current density cannot 
be uniquely defined. Consequently, for deriving fluid dy- 
namics and keeping the linear integration-operator struc- 
ture of the Boltzmann equation, the Tsallis formula is the 
one which provides us a straightforward derivation. 

Let us introduce the single particle phase-space dis- 
tribution function f k = f (t, x, k°, k) . The distribution 
function is normalized to the number of particles in the 



system, hence in an invariant phase-space volume element 
(x, x + dx, k, k + cik) at time t, the number of particles is 
given by, 

dhf (t) = f k d 3 xd 3 k, (9) 

where x = |x| and k = |k|. From now on we work with 
natural units and set k B = H = c = 1. 

The q-generalized entropy from Eq. ((TJ) can also be 
expressed with the help of the single particle distribu- 
tion function. This leads to the following relativistic q- 
generalized entropy four-current, 



S" = - dK fc"[(/ fc )*ln, (f k ) - f k ] 



(10) 



Here dK — gd 3 k/ ((27r) 3 fc°) is the invariant momentum- 
space volume with g denoting the number of internal de- 
grees of freedom, such as the spin degeneracy. The timelike 
component of the relativistic entropy four-current S — 
- / d 3 k (.g/(2vr) 3 ) [(f k ) q ln q (f k ) - f k ] retains the familiar 
Tsallis entropy, while the remaining components define a 
three vector. 

It is imperative to point out that in kinetic theory there 
are different definitions of the q-entropy four-current in 
use. These definitions are not equivalent with each other, 
and also lead to different thermodynamical relations as 
will be discussed later in this paper. For example, Lima 
and Silva et al. [2SU3D] applied S% = - f ^ Wfi \n q f k , 

while Lavagno [SJEl] used S£ = - J ^ Wfl [In, f k - 1]. 
These alternative definitions did not explicitly denote the 
finite size of the elementary phase-space cell volume (2irK) 3 . 
The corrections for quantum statistics can be added, both 
for fermions (z = 1) and for bosons (z = —1). A correct q- 
generalized quantum statistical entropy was proposed by 
Cleymans and Worku recently [33l[34ll35] . 



ft^ q fk + -{l-zf k ) q \n q {l~zf k ) 



(11) 



The z — limit corresponds to Eq. (TTU1) and was used 
in this form by Osada and Wilk g g H M gj]. 

The second law of thermodynamics demands positive 
entropy production, 



dK\a q {f k ) [fc^(/ fe ) 9 ] >0, 



(12) 



This entropy production formula suggests the following 
kinetic equation, first introduced by Lavagno 31 , from 
now on referred to as NEBE, 



k"0j k = C [f] 



(13) 



where C [/] is the collision integral and we use the nota- 
tion, 

fk = ifk) q ■ (14) 

This NEBE is postulated with the purpose to study clas- 
sical statistical effects due to correlations possibly induced 
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by finite system size, compared to the characteristic inter- 
action range. 

We also stress here that this distribution determines a 
different number of particles from Eq. ©, 



dM(t) = f k d 3 xd 3 k . 



(15) 



Of course this should not mean that the total number of 
particles in a system would turn out to be different after 
integration over the whole phase-space. 

The collision integral, C [/] , specifies the change in f k 
due to binary collisions among particles with initial incom- 
ing momenta of k'^ and outgoing final momenta p M , 
p lfl . It is constructed with the help of the invariant tran- 
sition rate W kk ,^. pp /. To fulfill the detailed balance prop- 
erty, the transition rate has to be symmetric with respect 
to the sequence of final states Wkkr^ppi — Wkki->prp, as 
well as symmetric for time reversed processes Wkki^ppi = 
Wppi^kki- The collision integral is formally written as, 

c \f\ = \j dK'dPdP'Wkk^pp, 

x(H g [f P J P >]-H q [f k ,f k ,]), (16) 

where the ^-generalized version of the assumption of molec- 
ular chaos, the so-called g-generalized Stosszahlansatz, is 
given as in [3"Tll3"2"] 

H q [f k ,fk>] = exp, [In, (/ fe ) + In, (f k , )] . (17) 

This formula is a postulate for the correlation between 
two particles with different momenta before and after the 
collision. For q = 1 it returns the familiar collision integral 
where the number of binary collisions around space-time 
coordinates is proportional to H\ [fk, fk'] = fkfk'- 

Clearly also other postulates can be made, which all 
satisfy an H-theorem [7]. The H q [fk, fk'} in the above for- 
mula, in general can be replaced by any fcfc'-symmetric 
functional. For example, using 



H q [f k ,fk'\ = exp g [In, (f k ) + In, (f k >) 
+ (l-g)ln, (f k ) In, (f k ,)] , 



(18) 



turns out to be the original Boltzmann type of factor- 
ization ansatz neglecting correlations, H q [fk, fk'] = fkfk 1 
even for q ^ 1. This however does not lead automatically 
to the classical kinetic theory, just then, if the k* 1 four- 
momenta are treated additively in the collisions. In this 
paper we assume that this is the case. 

Note that in the literature not only for the entropy 
but also for the NEBE collision kernel, slightly different 
versions of the kinetic equation and correlation function 
exist, involving the Ln, [27] or In,. [291130] functions. We 
return to discuss these later in the Appendix IA.11 

The entropy production vanishes in local g-equilibrium 
that is, = - J dK ln q (f ok ) C [/ ] = 0. Therefore the 

distribution function which satisfies c? M S" M [fo k ] = locally, 
is the ^-equilibrium distribution. This follows assuming 
that In, (/ofc) = a - fi^k^ hence, 



fok = exp (a - Pok^Ufj) . 



(19) 



Here ao and /3q = Pqu^ are collision invariants with u 11 
being a four- vector normalized to one, u^u^ = 1. These 
quantities are identified with the inverse temperature /3q = 
1/T, the chemical potential over temperature, ao = /x/T, 
and the fluid dynamical flow of matter u^. The formula in 
Eq. (fH)]) reduces to the well known Juttner [361137] distri- 
bution or relativistic Maxwell-Boltzmann distribution for 
q = 1, that is f Jk = exp (a - fiok^u^). 

We introduce the q-modified particle four-current and 
symmetric energy-momentum tensor, 



N» = dK Wfk, 



dK Wk v f k . 



(20) 
(21) 



Now, using Eq. (TT31 the four-divergences of the above 
quantities lead to, 



d^N" = J dK k^d^fk = J dKC [f] , (22) 
d„T^ = [dK k v k^dj k = [ dKk v C [f] . (23) 



Specifying the so-called collision or summation invariants, 
ip — a — f3^k^, the conservation of particle four-current 
and energy-momentum tensor are implied due to the con- 
servation of particle number and energy-momentum in in- 
dividual collisions. These follow from J dKaC [f] = and 
J dK^k^C [f] = for any solution of the Boltzmann 
equation, 



= 0, 



0. 



(24) 
(25) 



Having all conserved quantities at hand it is straightfor- 
ward to prove that the fundamental relation of equilibrium 
thermodynamics is given as, 



S" 



(26) 



where po is the equilibrium pressure. Projecting the above 
equation with the local flow velocity, we get the familiar 
form of the thermodynamical relation, 



s = -a n + A) (eo + Po) ■ 



(27) 



Here no = i7,(i.o) is the particle density, eo = J q (2,o) is the 
energy density and po = — ,7,(2.1) is related to the particle 

density as |^ = We introduced the following q- 



generalized thermodynamic integrals [25 : 
1 



J, 



K 



(2j + l)!! 
1 

(2j + l)!! 

q 

(2j + 1)!! 



dK(E k y~ 2j (A^k^KY fok, (28) 

dK (E k y- 2J (A^k^KY (fok) q , 

(29) 

dK (E k y~ 2j (A^k^KY (fokf 9 ^ 1 ■ 

(30) 
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where i,j > are natural numbers and the double fac- 
torial is (2j + l)H = (2j + 1)!/ (j\2 j ). Furthermore, we 
introduced = g^ v — u^u v , to project arbitrary four- 
vectors into other four- vectors orthogonal to v?. The par- 
ticle four- momentum is decomposed as k 11 = E\,u^ + k^ , 
where Ef, — k^u^ is the Local Rest Frame (LRF) energy 
of the particle and k^ = A^ v k v the LRF momentum. 
The LRF is defined by < flF = (1, 0, 0, 0). 

It is important to realize that the previously intro- 
duced thermodynamical integrals are consistent with ba- 
sic relations of classical thermodynamics. At fixed tem- 
perature (f ok ) q = (df Qk /da ) | ft, as well as q(fokf q ^ 1 = 
(d (fok) q /da ) \p 0) so one obtains 



fdl q( 



\ da J Pa 



K 



f dJq(i,j) 

V da 



(31) 



One shows by partial integration that the following recur- 
sive relations hold in equilibrium: 

1 i — 2j 

Jq(i,3) = ~~o- I q{i-hj-l) H a — ^(i-U)' ( 32 ) 

Po Po 

Kq(i,j) = ~!T^q(.i-l,j-l) ^ a^^i-hi)- ( 33 ) 

Po Po 



Furthermore, 



■Eq(i,j) — J^q(t,j)^0 — Jq{i+l,j)Pai (34) 
^q(i,j) = K-q(i,j)®a ~ fcq(i+l,j)Po- (35) 

These equations provide the ^-generalized Gibbs-Duhem 
relations. They not only reduce to the standard BG ther- 
modynamics for q — 1 but at the same time obey relations 
consistent with standard thermodynamics [251135] . Similar 
recursion relations have been obtained in framework of 
higher order generalized thermodynamics [38 , 39 , 40 , 41 . 



3 Out of q-equilibrium 

So far we have obtained a local q-equilibrium distribution, 
but yet we are interested in a more general solution of the 
NEBE. Assume that 



fk — fok + Sfk, 



(36) 



where Sfk denotes the deviation from local ^-equilibrium. 
It is customary to surmise Sfk <C fok ■ One should not mix 
this with an approximation (q — 1) -C 1, as pointed out 
in Rcfs. 22.23,24 . This distinction becomes physically 
important when considering the entropy production rate. 
Namely, the small (q — 1) approximation can spoil the 
positivity of the exact expressions in its leading order. 
This problem is based on the approximate behaviour of 
the basic deformed exponential: While 



e q (x) = cxp 



1 



In [l + (l-q)x] } >0 



(37) 



for any real argument, its small (q — 1) approximation, 

™2 



e q (x) 



l-(l-q) 



O 



(9-1)' 



(38) 



will have a negative sign for |af| > Furthermore, 

even when q > 1 would be postulated, in some formulas 
q* — 2 — q < 1 occurs - as it will be apparent later. There- 
fore such approximations loose the positivity property and 
should not be used in describing dissipative phenomena. 

The formula (|36[) can be approximated by a series 
expansion around local q-equilibrium. As previously dis- 
cussed, in local q-equilibrium we may define the collision 
invariant ip ok = ln q (f ok ) = ato — pWu M > hence 



fk (-0) = fk (tpok) 



dfk (Vtofc) 
dipok 



k + 0[<p 2 k 



hk + {hk) q 4>k + oU>i \ , 



(39) 



where 4> k — (ip — V'ofe) denotes the difference between the 
collision invariants. Therefore comparing the above result 
with Eq. (1361) we have to apply 



Sfk = fvk<t>k- 



(40) 



Because the primary physical quantities are associated 
with (fok) q henceforth Sfk has to be calculated too. To 
do this, we expand f k = fok + Sfk around f okl 



fk = (fok + fak<Pk + O [<j)i_ 



= hk+q{hkf q ' 1 <f>k + 0[^ 



and so we obtain, 



Sfk = q (fok) 2q 1 <j>k 



(41) 



(42) 



Let us recall the q-generalized Stosszahlansatz ([TTf . To 
further simplify our discussion we expand it in series, 



Hq [fk, fk 



Hq [fok, fok'] H 
dH q [/ofe, fok'] 



of\ 



Ok' 



dH q [/pfc, f ok ' 
dfok 

Sfk' + O [Sf 2 ] 



Sfk 



where 



dH q [f k, fok' 



Of 



OA- 



dHq [f 0k , f 0k ' 



d.f\ 



OA:' 



— (Hq [fok, fok']) q (/ofc) 
= (Hq [fok, fok']) q (fok>) 



(43) 

(44) 
(45) 



Now, making use of Eq. (1401) we arrive at, 

Hq [fk, fk') = H q [f 0k , fok>] + (H q [f 0k , /ofc']) 9 (4>k + <f>k') ■ 

(46) 

In this way the collision integral from Eq. (|16|) is approx- 
imated up to first order in deviations from q-cquilibrium 

as, 

C[f] = C [f ] + C [Sf] . (47) 
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Here C [/ ] oc H q [f 0p , f 0p '} - H q [/ofe, /ofe'] vanishes by def- 
inition, H q [/op, /o P '] = H q [fokJok'}. Then C[f] = C [Sf] 
is given by the integral, 

C[5f] = \\ dK'dPdP'W kk ,^ pp , (H q [f Qk Jok']) q 

x {(j) p + 4> p > - 4>k - <j>k'), (48) 

and so the NEBE (O simplifies to, k»dj k = C [Sf]. 

Using the out-of q-equilibrium distribution function, 
the particle four-current and the energy-momentum tensor 
from Eqs. (|20I21[) take the general form, 



matching conditions, 



5n = J dKE k 5f k = 0, 
Se= JdKE 2 k Sf k = Q. 



(59) 
(60) 



The isotropic pressure p separates into two parts: p 
Po + II where the bulk viscous pressure becomes, 



n 



1 



dK{A a ^k a kp)5~f k . 



(61) 



A" 1 = nil** + V", 



(49) 



T» v = eu^u v - pA"" + u"W u + u v W^ + tt"" . (50) 

We can uniquely identify the fundamental fluid dynami- 
cal quantities like the particle density, energy density and 
isotropic pressure by following the matching procedure, 



n = u M JV" = J dKE k f k , (51) 

e = Utl u v T^ = JdKE 2 J~ k , (52) 

P = -\A„ V T^ = -±JdK (A a ^k a kp) f k . (53) 

Similarly the particle diffusion and energy-momentum four- 
currents as well as the stress tensor arc defined as, 



The conservation equations for Nq and Tq U turn to the 
Euler equations for a perfect g-fluid, while Eqs. (|24I25I) 
together with Eqs. (|49l IBTJl) constitute the equations of 
dissipative g-fluid dynamics. 

Although we motivated the g-fluid dynamical equa- 
tions form the NEBE, these constituting equations can 
also be postulated as the g-generalized version of the con- 
servation equations based on the non-extensive g-entropy 
four-current and the corresponding laws of thermodynam- 
ics. The equations of perfect g-fluid dynamics are closed 
by an Equation of State (EoS). Without a first order phase 
transition it implies a vanishing local entropy production 
in g-equilibrium, O^Sq = 0. However, the equations of q- 
fluid dynamics presented so far in this paper are not closed 
because the out-of g-equilibrium f k is unknown. Therefore 
in this case one needs additional relations beyond the EoS 
for completing the system of equations. 



V 1 ' 



A^N C 



dKk^Sf k , 



(54) 



4 The Navier-Stokes equations and transport 

W = A» a u p T a P = / dKE k k^5~f k , (55) coefficients 



ntt u = T <j*u) = / dKk^k v) 8f k 



(56) 



where Sf k = f k — fok and we used the following con- 
straint / dK (E k ) 1 (A^kf.kv)" 1 k^ 1 ...fc^/ofe = for 
any l,n,m natural numbers. Here T^") = A' lva ^T a p, 
with A^ vafi = \ [AP ,ol A^ v + A'^A 13 ^) - \A^ ' A afi pro- 
jector, selecting out the traceless, symmetric and orthog- 
onal to tin part of T^". 

The particle four-current and the energy-momentum 
tensor calculated from Eqs. (|20I21I) by substituting the 
local ^-equilibrium distribution function from Eq. (|19p . 
leads to the so-called perfect q-fluid decomposition, 



K 



(57) 
(58) 



Taking advantage of the previously derived approxima- 
tions we reduce the problem to finding Sf k = fo k (f> k near 
local g-equilibrium. There are several well known approx- 
imations to obtain the deviations from local equilibrium 
as well as to calculate the transport coefficients. Here we 
follow the traditional approach of Chapman and Enskog 
[101411142) and apply this method for the NEBE. To this 
end we note that this method is iterative, and here we 
take only the first step, which results in the well known 
Navier-Stokes equations of fluid dynamics. 

Assuming that in the vicinity of local ^-equilibrium 
Sfk <C /ofc, the non-equilibrium contributions from the 
streaming term are also assumed to vanish k^d^5f k = 0. 
Therefore k^d p f k ~ k^d p f ok , i.e., 



where no = 
e = u^uJT^ 



J dKE k fo k is the particle density, 
= / dKElf ok is the energy density and 

Po = -§A^ T cT = -lfdK (A<*Pk a kp) f ok is the pres- 
sure in ^-equilibrium. The particle density and energy den- 
sity are usually assumed to be unchanged from their equi- 
librium values n = uq and e — eo, also cited as the Landau 



P o ~(k a k A af3 )+E k a o -E 2 k j3 o 



k^djok = qfok 1 

+ (1 - E k h^) fc<">V M ao - fok^k^o-^ 



(62) 



Here, 6 = V M u M is the expansion scalar, a^ v — V^?W 
is the shear-stress tensor and V^ap is the gradient of 
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the chemical potential over temperature. The proper time 
derivatives do and $q can be expressed as 

d = "of ~(2,0) (^0^(2,0) - £g(3,0)) (63) 
$0 = ~(2, ) (^o£ g (i,o) - £z(2,0)) O, (64) 

where ho = (e + Po) /n = £9(3,1) /£g(2,i) is the enthalpy 
per particle and T> q{lJ) = K, q{l _ 1J} IC q{l+lj) - £g (iJ) . 

The main equation we need to solve in the first order 
Chapman-Enskog method reduces to, 



k»dJ k = C[5f} 



(65) 



where the left- hand-side is given by Eq. (R)2"1) . while C [Sf] 
is given by Eq. (|48|) . Inspecting this integro-differential 
equation, we search for its solution using the ansatz |41j . 

<Pk = (<Po + ViE v + ip 2 El) 6 + (<p a + (p 4 E k ) k^V ^ 
+ ^ 5 k^k^a^, (66) 

where the ipiS are yet unknown coefhcients. In an earlier 
work [25] we calculated the transport coefficients based 
on the Anderson- Witting model [43] where <pk was given 
directly by Eq. (j52"j) using the so-called relaxation time 

approximation, C [Sf] = -E k (f k - f ok j /tc, with t c be- 
ing the mean time between collisions. Similar calculations 
were also done in the non-relativistic limit by Bezerra et 
al. [2] using an alternative NEBE from Refs. [251131]] . 

Here we proceed with the calculation of the (fi coefh- 
cients. Using Eqs. (|59I60[) and the definition of the bulk 
viscous pressure from Eq. (j61"j) together with Eq. (|42|) we 
get three coupled equations for the first three coefficients: 



Sn = 



q J dKE k (fok) 2 ' 9k 

= (<A)£g(l,0) + ^1^9(2,0) + ^2^(3,0)) 

5e = qJ dKE 2 k (/ofe) 2 ^ 1 A 

= (<A)£g(2,0) + V?l£g(3,0) + V2^3(4,0)) ^> 

n = 



\2g-l 



-| J dK{A a f > k a k [j ){f ok f q ~ 1 9k 

' (<A)£g(2,l) + £<i(3,l) + ¥>2£g(4,l)J 



(67) 
(68) 
(69) 



For Sn — Se — from the matching conditions (|59I60[) . 



one concludes that 



<P2 



D n 

n (£. 



9(2,0) £g(3,0) 



£9(l,0)£g(4,0) ) 



D 



11 



nv, 



9(2,0) 



6 D n ' 
where the denominator is given by 

Dn = £9(2,1)^9(3,0) + £g(4,l)£ ) g(2,0) 

+ £9(3,1) (£9(2,0) £9(3,0) — £g(l,0)£g(4,0)J 



(70) 
(71) 
(72) 

(73) 



It is important to note that the matching conditions are 
not unique and other choices may be made [45ll46ll47] . 
In general one may assume that Sn ^ Se ^ 0, so the 
coefficients (fi (6n,Se,II) are functions of all scalar non- 
equilibrium corrections. However, as we shall see it later, 
the constraints to the NEBE only supplies (f2 , the other 
two coefficients are left undetermined. Therefore, in gen- 
eral we need additional input to specify the non-equilibrium 
scalar corrections, and in this paper we work with the 
matching conditions from Eqs. (|59I60[) . 

Similarly, the definitions of the particle and energy- 
momentum diffusion four-currents and the stress tensor 
lead to, 



V = qj dKkM (fokf^ 1 <}>k 

= (V3£g(2,l) + ^£9(3, 1)) V M a , (74) 
W» = qJ dKE h kV> {hkfi- 1 (f>k 

= (<^3£g(3,l) + ¥>4£g(4,l)) V M a , (75) 

7r^ = qJ dKk^k^ (/ 0fe ) 29 "Vfc 

= 2< y35 £9(4, 2 )^ Il/ . (76) 

Here we used, k^k^ = k^k^ - \ {A af3 k a k ) A* v , and 
the general orthogonality relation, 



S( 1/1 ...k Vn } 



J dK F(75 k ) k^ 1 ...k^k { 
7^%T^::r / dKF(E^) (A^k a k ) m , (77) 



with F(£'k) being a function of energy |40) . 

These results for the dissipative quantities present the 
relativistic Navier-Stokes relations, usually written in the 
more familiar form, 



n = -ce, 



K 



'T 



hoP 2 
; " - 2ria"". 



Va , 



(78) 
(79) 

(80) 



The relativistic counterpart of the Stokes result, 77 = — (9, 
relates the bulk viscous pressure to the expansion rate 
by introducing the coefficient of bulk viscosity £. The 
Fourier-Navier- Stokes law in Eq. ([79]) relates the heat-flow 
= — hoV^ to the temperature gradients in the sys- 
tem, V^a = -h T- 2 (V^T - TO 1 ) with the coefficient 
of heat- or thermal-conductivity k. Finally the Newton- 
Navier-Stokes relation between the stress and shear is given 
in Eq. (|80p with the coefficient of shear viscosity n. 

The tfi coefficients are obtained inserting <p k from Eq. 
(f6"6"]l into the collision integral (|4"8"|) and matching the gra- 
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dients on both sides. The collision integral leads to, Similarly, the coefficient of heat-conductivity can be 

r- calculated, by matching the gradient V^ao on both sides 

C [Sf] = - dK'dPdP'Wkkr^ppi {H q [f ok , fok'}) q of the integrated NEBE. This leads to, 

x[^(E 2 p +E 2 pl -El -El) 9 _ /C g(3 ,D - V ^.p 

— v • l 9i ; 

As 



+^ A (E p p^ + E p ,p'^ - E k k^ - E k ,k'^) V^ao 
+(p 5 (p ifi p^ +p'^ P ' u) -k^k v) -k'^k' u A o pv 



The </?3 coefficient can be found only after fixing the LR 
frame. Here we have at least two choices. We may choose 
(81) V 1 = 0, which corresponds to Eckart's definition [38] . 
and hence the fluid dynamical flow velocity is fixed to the 
where the contributions vanish identically for ip , Vl and flow of con served particles, = A^/y^V Another 
<p 3 due to the conservation of charge, energy and momenta choice due to Landau and Lifshitz g5] defines the flow as 
in binary collisions. ^he time-like eigenvector of the energy-momentum tensor, 

For later purposes we introduce the tensor, ^ = T ^uJ JT^u^T.^ . This choice is equivalent to 

symmetric upon the interchange of indexes (u, u) as well w „ = Q choosing the Landau framej ne arrives at 
as (a, p), while traceless for the latter, X^ ' g a p = 0: 

1 r _ ^-9(4,1) / no x 

= ~J dKdK'dPdP'W kkl ^ ppl (H q [f ok , f ok ']) q ^ - '^K^' (92) 

x WW (p Q / + p' a p ,f) -Wk? - k' a k ,fs ) . (82) an d thus the coefficient of heat conductivity becomes 

Furthermore, one can show that X^ al3 is generally de- D - 3 ^ 

composed as, K — fi^oPoip . (93) 

<S(3,1) 

X^ a ' 3 = (XmW + X 2 A^) (u a u? - ^A a A Finally, the coefficient of bulk viscosity can be ex- 

^ ' pressed with the help of Eqs. (|S3llM|) : 

+ AX 3 u^A v ^ a u^ + X A A^, (83) 

p-i 

with the coefficients, ^ = _ [p ]C q{i<1) V q{m - n V q(m 

X l =X^ a ^u u u a up = -X^ a ^u v A a p, (84) u (r r r r W «xa\ 

1 1 -no/10 ^,(2,0)^(3,0) _ ^-9(1,0)^(4,0) J J ) (94) 

A 2 = ~X» ua A^u a u p = --A^ Q <VA*/3, (85) hence from Eqg ^ ^ we obtain ^ 

X 3 = \x^ U{lJl A v){a u &) = \x^u^A va u^ (86) DjJ 

? d C = ^2^ ' ^ 

X 4 = ix^^. (87) _^ 9(2 ^ 0) 

These results formally resemble their classical coun- 
Applying these results and integrating Eq. ([55) we get terparts within BG statistics. An explicit g-dependence 
three equations which are proportional to the different occurs in the thermodynamic integrals only while the co- 
Navier-Stokes type gradients occurring in Eqs. J78J80|). efficients X 1 ,X 2 ,X 3 and A 4 differ in their arguments. 
For example, collecting terms proportional to a a/3 leads i n the fallowings we attempt to simplify the collisional 
to, integral. Recall the H q [fk, fk'] functional from Eq. (flTl) 

and calculate it in equilibrium getting 



-f3 j dKWWk (a k p) (qflr 1 ) 
^ J dKdK'dPdP'W kk ,^ pp , (H q [/ 0fe , /, 



n!.'\) q 



(H q [f ok , fok'}) q = (cx Pg [In, (/ofe) + hi, (f ok >)}) Q 

= [l + (l-q) (yjok + ifov)]^ ■ (96) 



xWk 



v (p{uPp) +p\ a P'i3) _ k(akp) — ^{a^'f})) > (88) Now ma king use of Eq. ^ one concludes that 



which briefly reads as, -2/3 /C g(4: 2)^ Q/3 = ^X^^ . H g Ifok, fok>] = exp g {ip Qk + ip ok >) . (97) 

From this we easily get that . 

lhe above formula for g = 1 leads to the classical re- 

_ _ 2 fl ^( 4 - 2 ) / 8q \ suit, ffi[/ofe,/ofe'] = fjkfjk' with / Jfc ee exp(i/>ofc) = 

^ 5 X A ' expfctn - BnWu,.). Substituting this result into Eq. (1521 

t. i • ,1 • , • *T??5r i ... we arrive at, 

Replacing this result into Eq. (|76p and comparing it to 

Eq . (1501) leads to the coefficient of shear viscosity, . 1 r 

2 X^ al3 = -J dKdK'dPdP'Wkk^pp, (ex Pg (Vofe + V'ofeO) 

V= <P5fc q (±,2) = -2fi Q -^-. (90) xWk" (p a pP +p' a p' 13 -k a kP -k' a k'P) . (98) 



8 



T. S. Biro, E. Molnar: Non-extensive statistics, relativistic kinetic theory and fluid dynamics 



This integral can be rewritten with the help of a Man- 



delstam invariant s = (fc M + k'^) = (p^ 
transition rate is given as |40j . 



P 



So the 



W kkl ^ pp , = (2tt) 6 sa (s, 9cm) 5 i {k» + k'» - - p" 4 ) . 

(99) 

Here a (s, 8cm) is the differential cross section, 9cm is the 
scattering angle in the center of momentum frame, while 
the J-function represents the energy-momentum conserva- 
tion in binary collisions. 

Substituting the transition rate into Eq. (j98|) . we get 
the following expression in the center of mass frame for 
the p-dependent integral: 



l p (a T ,s,m) 



1 



d 3 p 

(P°f 



(100) 



Here we defined the total cross-section <tt (s) as the inte- 
gral over the differential cross-section and the solid angle, 
dQ = dip J* sin OcudOc m, hence 



i r 

a T (s) = - J dQ a (s, 9 C m) ■ 
Finally the X^ a P tensor can be given as, 



(101) 



X v, a p = - i dKdK , Tp ((TT; Sj m) ( exp ^ ^ m + 



x k^k v l - 



pap/3 _ 



Am 



a/3 



pa p/3 



(k a k P + k' a k'P) 



(102) 



where we introduced the total momentum P M = k^+k'^ = 
pfj- +p'v. This integral, contrary to the classical examples 
is not easy to handle analytically due to the fact that 
only exp g [x + y + (1 - q) xy] = exp q (x) exp q (y), hence 
the k integrals can not be factorized into thermodynamical 
integrals. 

As already shown the thermodynamical quantities and 
transport coefficients change with g-parameter. The fact 
that the collision integral does not factorize, signals a di- 
rect g-dependence and the effect of the correlations in- 
cluded in the NEBE. Therefore it also follows that the 
original Boltzmann equation {q = 1) neglects any kind of 
memory effects and two particle correlations, hence the 
transport coefficients calculated from such an equation 
also lack correlations. 



5 Conclusions and outlook 

In this work we derived the ^-generalized versions of the 
classical Navier-Stokes equations of relativistic dissipa- 
tive fluid dynamics from a g-generalized Boltzmann trans- 
port equation. These equations were found based on the 
Chapman-Enskog method. 



We showed that starting from a (/-generalized trans- 
port equation, it is justified to apply standard methods to 
calculate transport coefficients. These calculations lead to 
relations for all transport coefficients formally similar to 
those one would obtain using the traditional Boltzmann- 
Gibbs distributions. The main difference is contained in 
the recursive rules for the thermodynamic integrals. 

However, unlike in traditional fluid dynamics, the ten- 
sorial collision kernel X^ va P does not factorize into prod- 
uct of simple thermodynamical integrals. This remarkable 
property is a consequence of the g-deformed exponential 
function describing the extended local g-equilibrium dis- 
tribution. 

We also add that the derivation presented in this work 
is applicable to slightly different kinetic equations as dis- 
cussed in Appendix IA.1I as well as to the one suggested 
by G. Kaniadakis [3,50,51 . Furthermore, we also show in 
Appendix lA.2l that one obtains causal fluid dynamics from 
the NEBE by using Grad's method moments. 

These methods applied to the ^-generalized Boltzmann 
equation extend the applicability of dissipative fluid dy- 
namics for q ^ 1 by including long range interactions and 
correlations, but for example the resulting Navier-Stokes 
equations are still parabolic, hence problems related to 
acausality are not solved by introducing a non-extensivity 
parameter q in this framework. Therefore this means that 
the causality problem is not rooted in the non-extensivity 
of entropy, but the entropy should be extended to include 
dissipative quantities which is well known from irreversible 
thermodynamic theories [52,53 . 

On the other hand, we remark that Osada and Wilk 
inRcfs. [22,23,24 associated the perfect g-hydrodynamics 
with the classical Navier-Stokes equations using a so-called 
non-extensive/dissipative correspondence (NexDC). This 
approach makes a direct correspondence between the purely 
g-dependent conserved quantities, Eqs. ([571 1551) for q^l, 
and the general dissipative structure as given in Eqs. (|49l 
I50p for q = 1. Therefore the dissipative quantities are iden- 
tified with differences given by perfect g-fluid dynamical 
quantities for q =^ 1 and classical dissipative fluid dynam- 
ics for q = 1. 

This novel method may reveal the correlations induced 
by long range effects that in general contribute to dissipa- 
tion and entropy production. However our approach goes 
beyond by generalizing the fluid dynamical equations to 
out of g-equilibrium states, hence in principle one can di- 
rectly estimate the g-parameter without conjecturing a 
NexDC. 
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A Appendix 

A.l Remarks on a different NEBE 

Here we briefly discuss different NEBEs advocated by 
some authors. The existence of different kinetic equations 
shows that there is an ambiguity about the correct equa- 
tion of motion, although they are consistent with the same 
g-generalized entropy. 

Let us recall the entropy four-current from Eq. (TR 
and rewrite it in equivalent forms: 



dK fcq/ fe hv (A)-/*], 
dK k" [f k Ln q (f k ) - qf k ] . (103) 



S" 



<1 

The requirement of positive entropy production turns out 
to be, 



and 



q / dKln q , (/ fc ) [fc^/fc] > 0, (104) 



dKLn q (f k ) [k^fk] > 0. (105) 



The above formulas suggest a different version of the q- 
generalized Boltzmann equation. This alternative NEBE 
can be written as, (denoted with hats to avoid confusion) 



Wdufk =C f 



(106) 



Here the collision integral is 
1 



C 



dK'dPdP'W kk ,^ pp , 



x H, 



fp i fp' 



H„ 



fk i fk' 



and 



fk j fk' 



exp 



Ln g ( f k ) + Ln 9 ( f k , 



(107) 



(108) 



The outcome of the H-theorem as shown by Abe [37] , is 
equivalent to d^S^ = — J dKLn q (fok) C [/] = 0, whence 

the collision invariant is tfjQ k = Ln 9 (^fokj , leading to 



f ok = E q (a - PqWu^) = exp ? , 



(109) 

We further note that Lima and Silva et al. [2j), ; .3j]] used 
a slightly different g-generalized Stosszahlansatz given as 
H* q [ft, ft] = exp, [hv (ft) + ln r (ft,)}. Actually, this 
can be achieved by the duality transformation, ln g * (x) — 
^Ln g (x), and so their equilibrium distribution function 
becomes, 



(110) 



It is clear that the different NEBE's lead to different sta- 
tionary solutions. However the cited solution of Abe as 
well as of Lima and Silva et al., consequently leads to dif- 
ferent thermodynamical relations than ours. 

All these different ^-equilibrium distributions can be 
used to define different particle four-currents and sym- 
metric energy-momentum tensors. For example, denote 

W (fok) = K, (; ofe ) = fff" and (/ ofe ) - 

whence by using Eq. (|103|) together with the correspond- 
ing stationary solution (|109[) . we obtain 



st; 



Oo 

q 



dKk^fok + — I dKk»E k f ok + / dKWf 



q 



I Ok 



q q 



(111) 



Projecting the above equation and introducing &q — ao/q 
and /3 = po/9 we get, s = -a h Q + /3 e + h . Cor- 
respondingly one has to define, the following thermody- 
namic integrals, similar to Eqs. (I28I29I30[) , 



(112) 



J, 



, 2/ ._ , dK(E k y- 2j (A'- k r k r )'l\ lh . 



K 



(2j + 1) 



dK(E k y- 2j (A^k^Kf (/< 



>0k 



(113) 

2-q 

(114) 



Now we rewrite the thermodynamic relation as, 

Sq = —&0 Jq(l,0) + PoJq{2fl) + Jq{l,Q)- (H5) 



However at fixed 

2-q 



temperature f ok = (d (fok) /<9a ) | /3o 



and q [fok] = (dfok/dao) \p hence 



J, 
while 



dl q{i>j) \ ~ Id J, 

, K, 



, (116) 



h(^3) ~ _ o-4(i-lJ-l) + a ( 117 ) 
Po PO 

These relations actually lead to an equilibrium pressure 
given as, po = —J q (2,i) — Po~ ^9(1,0) an d so the last term 
from Eq. (|115[) cannot be identified with the pressure 
times inverse temperature $oPo Jq(i,o)- Therefore the 
equilibrium state defined by Eq. (|109[) and Eq. (jllOp . are 
not consistent with the ideal gas EOS in its classical form, 
while the one utilized in Eq. (fH?]) is, (3 po = Zq(i, o)- T his 
docs not necessarily mean that the NEBE in Eq. (1106[) or 
its stationary solutions are ruled out, but one has to keep 
in mind that they lead to weird thermodynamic relations. 
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Similarly, as previously presented we can define the 
deviations from equilibrium as, 



1 / « \2- 9 A 

fk = fok + - yfok j 4*k + 



where 4> k = - ■tpok) and so Sf k = \ {jok\ 4>k- The 
collision integral also simplifies and formally corresponds 
to Eq. flUJl, 



of pr hence also of H, q^ 1 , and tt^ . The relaxation 

of the dissipative quantities towards their Navier-Stokes 
values is given with time scales given by the corresponding 
relaxation times Tn, T q and t w . The calculation of the 
collision integral actually involves the very same procedure 
as presented in the previous sections, while the relaxation 
equations are also given in exactly the same form as in 
the classical case. Once again the 1 modifications are 
embedded in the thermodynamical and collision integrals. 



C[Sf] = \ I ,1K',W,W'\Yi L , ,,„ 



fok, JOk' 



U> P + 4> P ' - 4>k — 4>k'^j ■ 



(118) 



Therefore independently of the underlying NEBE the 
transport coefficients can be calculated by any classical 
method of choice. However as we have shown, differences 
arise in the thermodynamical integrals and the EoS. 



A. 2 Grad's method of moments 

Here we discuss an alternative method to calculate the 
transport coefficients. This method, originally due to Grad 
[51] leads not only slightly different transport coefficients 
but also to different equations of motion for the dissipative 
quantities. This is due to the different choice of parame- 
ters for the non-equilibrium distribution function and the 
particularities of the method, see for example Refs. [38l 
HP1HT1I5HI551I5W571I55] ■ for more details. 

Here we follow Refs.|57j[55] and introduce the following 
irreducible tensor moment for the deviation form local q- 
equilibrium, 



rr 



dKElk^ 1 



.k^Sf k , (119) 



where the index £ indicates the rank of the tensor such 
that I = corresponds to the scalar p r , with power r of 
the energy . Similarly one can introduce irreducible mo- 
ments involving the solution of the other NEBE from Eq. 
(flTOl) . i.e., p'f 1 -^ = fdKElk^ ...k^Sfk. The gen- 
eralized irreducible moments from Eq. (|119[) are identified 
with the dissipative quantities via the following relations, 



5n = =0, Se = p (2 ) =0,11 



-P(0), 



7" 



p 



(())' 



-LIV 



(120) 



Now, rewriting the NEBE (TTH)) in the following form, 



dr dr 

E^k v V v Sf k +E^C[Sf] , (121) 



we obtain exact equations for -rrpi 1 * 1 This method is 



presented and analysed in great detail in Refs. [571158] . 

The main results of Grad's method are the so-called 
relaxation equations, which determine the time evolution 
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